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I The dominant theoretical uncertainty in extracting |V(d/V(s| from the ratio of branching 
^ \ ratios R = B{B {p,u)'y)/B{B K*'y) is given by the ratio of form factors ^ = 

re-examine ^ in the framework of QCD sum rules on the light- 
cone, taking into account hitherto neglected SU(3)-breaking effects. We find ^ = 1.17 ± 
0.09. Using QCD factorisation for the branching ratios, and the current experimental 
average for R quoted by HFAG, this translates into = 0.192 ± 0.014(th) ± 

0.016(exp). This result agrees, within errors, with that obtained from the Standard Model 
unitarity triangle, | Vi(i/Vis|g^^ = 0.216 ± 0.029, based on tree-level-only processes, and 
with \Vtd/Vts\Am = 0.2060t[!;[J[]^J(th) ± 0.0007(exp), from the CDF measurement of 
oscillations. 

This version differs from the original version of the paper, published as JHEP 
04 (2006) 046, by the inclusion of the new BaBar measurement oi B — > p{^)l 
presented at ICHEP 2006, which significantly shifts the results for |Vid/T4s|- 



*Patricia. Ball@durham.ac.uk 
tRonian.Zwicky@durham.ac.uk 



1 Introduction 



Recently, the Belle collaboration measured the b ^ d penguin-dominated decay B — > 
(p,a;)7 [1], whereas BaBar obtained an upper bound in 2004 [2] and presented a measure- 
ment at ICHEP 2006 [3]. Assuming the Standard Model (SM) to be vahd, this process 
offers the possibility to extract the CKM matrix element \Vtd\, in complementarity to the 
determination from Bd mixing and the SM unitarity triangle based on |V^„6/V"c6| and the 
angle 7. In order to extract \Vtd\ from the measured rate, one needs to know both short- 
distance weak and strong interaction effects and long-distance QCD effects. Whereas the 
former can, at least in principle, be calculated to any desired precision in the framework of 
effective field theories, and actually are currently known to next-to-leading order in QCD 
[4], the assessment of long-distance QCD effects is notoriously difficult. After a long history 
of phenomenologically or l/A^c-motivated factorisation formulas, QCD factorisation [5, 6] 
has provided a consistent framework allowing one to write the relevant hadronic matrix 
elements as 



Tf^^(0)7;^+ f diduTl\i,u)ct>B{i) 
Jo 



(1) 



Here e is the photon polarisation 4-vector, Qi is one of the operators in the effective 
Hamiltonian, T^-^^ is a S — > y transition form factor, and 0s, 0y;_L are leading- twist light- 
cone distribution amplitudes of the B meson and the vector meson V , respectively. These 
quantities are universal non-perturbative objects and describe the long-distance dynamics 
of the matrix elements, which is factorised from the perturbative short-distance interactions 
included in the hard-scattering kernels T/ and T/^. The above QCD factorisation formula 
is valid in the heavy-quark limit m^, 00 and is subject to corrections of order Kq^qd/itli). 
Although it is possible to determine \Vtd\ from the branching ratio of -B ^ {p,uj)'~f itself, 
the associated theoretical uncertainties get greatly reduced when one considers the ratio 
of branching ratios for B — > K*^ and B — > {p,uj)j instead. One then can extract | Vtd/V(s| 
from 



B{B^{p,u;h) 
B{B K*^) 



td 



V, 



ts 



ml./ml) \Tr(0)J l^ + ^^l- 

where the estimates of AR available in the literature lie between, approximately, and 
0.2 [5, 6]. AR contains all non-factorisable effects induced by T^'^^ in (1). As \Vts\ = |Kfe| 
in the SM, up to a small correction ~ 2%, and \Vcb\ is known with a precision of 2% [7], 
\Vtd\ follows immediately from |\4(i/Vts|. The theoretical uncertainty of this determination 
is governed by both the ratio of form factors T/^*(0)/Tf''^(0) and the value of AR, which 
parametrises not only SU(3)-breaking effects, but also power-suppressed corrections to the 
QCD factorisation formula. The aim of the present paper is to re-examine the size of 
SU(3)-breaking corrections to T/^*(0)/Tf''^(0) from QCD sum rules on the light-cone and 
to determine a value of |14d/T4s| from (2), evaluating AR in QCD factorisation; we will 
address the issue of power-suppressed corrections to AR in a separate publication [8] . 

Our paper is organised as follows: in Section 2 we discuss the QCD sum rule for the 
ratio of form factors T/^*(0)/Tf''^(0) and its dependence on SU(3)-breaking parameters. 
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In Section 3 we extract a value of |V^d/Vts| from the experimental branching ratio, using 
QCD factorisation for the calculation of AR. We summarise in Section 4. The appendix 
contains some formulas relevant to the calculation of NLO evolution of twist-2 vector- meson 
light-cone distribution amplitudes. 



2 The Form-factor Ratio Tf^^^/T^^^ 

In this section, we present a concise formula for the form factor Tf^^{0), as obtained 
from QCD sum rules on the light-cone, and discuss the hadronic quantities that enter this 
expression. We do not discuss the technique of QCD sum rules itself, or that of QCD sum 
rules on the light-cone, for which we refer to the literature [9]. Suffice it to say that the 
light-cone sum rule for Ti is based on the light-cone expansion of the correlation function 
of the chromomagnctic dipolc operator Qj and the interpolating field 5^756 of the B meson. 
The expansion is in terms of the convohition of process-specific pcrturbativc kernels and 
universal meson light-cone distribution amplitudes (DAs) of the final-state vector meson, 
which are ordered in terms of increasing twist. These DAs have been studied in Refs. [10, 
11], mostly for the p meson, including two- and three-particle Fock states up to twist 4. 
An extension to the K* meson is in preparation [12]. The light-cone expansion is matched 
to the description of the correlation function in terms of hadrons by analytic continuation 
into the physical regime and the application of a Borcl transformation, introducing the 
Borel parameter and exponentially suppressing contributions from higher-mass states. 
In order to extract the contribution of the B meson, one describes that of other hadron 
states by a continuum model, which introduces a second model parameter, the continuum 
threshold sq. The sum rule then yields the form factor in question, Ti, multiplied by the 
coupling of the B meson to its interpolating field, i.e. the B meson's leptonic decay constant 
/b- At tree level, the sum rule for ^^^"^^(O) then reads, to twist-4 accuracy: 



m 



rrib 



f 



2u 



J uo 



2u 2^^ ^ ' 8u 



1-u 



du 



Id r . r , u-ai , ,, . , 

;— / dai / da2 ^ ^ -^(^J + Via, 

uduJo Jo ^oti * 



uu(j)_i{u) + 2Il{u) + uH3{u) 



mi Juo 

- r dai r da2 — {S{a) - S{a) + ^ {a) - T^'^ (a) + T^^^ {a) - T^^ {a 

Jo Jo «3 ^ 



1 . ^ ^ 



(3) 



2 



nib 



Juo ^6 \ ^6 / 



(4) 



where Uo is given by ml/sQ. fy and fy are the decay constants of, respectively, longitu- 
dinally and transversely polarised vector mesons. $, g^^'"^\ II and are two-particle 
distribution amphtudes and integrals thereof, as defined in Ref. [13]. A, V, 5", S and t/^^ 
are three-particle DAs. u is the longitudinal momentum fraction of the quark in a two- 
particle Fock state of the final-state vector meson, whereas q;i,2,3, with ^cti = 1, are the 
momentum fractions of the partons in a three-particle state. The light-cone expansion is 
accurate up to terms of order (my/m^,)^. Although we only write down the tree-level ex- 
pression for the form factor, radiative corrections are known for Ri [14] and the two-particle 
contributions to R2 [13], and will be included in the numerical analysis. All scale-dependent 
quantities are calculated at the (infra-red) factorisation scale /^^ = m\ — ml- The form 
factor itself carries an ultra-violet scale dependence, which however cancels in the ratio. 

It is clearly visible from the above formula that the respective weight of various contri- 
butions is controlled by the parameter my /nib; the next term in the light-cone expansion 
contains twist-3, -4 and -5 DAs and is of order {mv/nibY- Nonetheless, (4) cannot be in- 
terpreted as l/nib expansion: for nib 00, the support of the integrals in u also becomes 
of 0{l/nib), as 1 — -Uq ~ 1 — ml/so ~ Uo/nib, with cjq ~ 1 GeV a hadronic quantity [15]. 
In this case, the scaling of the various terms in nib is controlled by the behaviour of the 
DAs near the end-point m — > 1. For finite nib, however, the sum rules are not sensitive to 
the details of the end-point behaviour, as we shall see below. Numerically, the expansion 
in terms of my / nib works very well and is a reformulation of the ordering of contributions 
in terms of the parameter 5 introduced in Ref. [13]. 

We have already discussed Ti in Refs. [14, 13]; in the present paper we focus on the 
ratio 

which governs the extraction of |Vts/Vtd| from B V'~f decays. Our sum rules can of 
course be used to determine each form factor separately, but we expect the ratio to be 
more accurate, because ^ is independent of the 5-meson decay constant Jb and also, to 
very good accuracy, of nib a-nd the sum rule parameters and Sq] we shall come back 
to that point below. Hence, in this paper, we will not re-analyse the absolute values 
nor, consequently, the branching ratios themselves. However, ^ is very 
sensitive to SU(3)-breaking effects in the DAs, and it is precisely these effects we shall 
focus on in this paper. A similar analysis for the ratio of the D ^ K and D ^ n form 
factors was carried out in Ref. [16]. 

Compared with our previous results of Refs. [14, 13], in this paper we implement the 
following improvements: 



• updated values of SU(3)-breaking in twist-2 parameters; 

• complete account of SU(3)-breaking in twist-3 and -4 DAs; 
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• estimate of higher-order conformal contributions to twist-4 DAs, using the renor- 
malon model of Ref. [17]; 



• NLO evolution for twist-2 parameters. 

Before presenting numerical results for ^, let us first discuss the values of the hadronic 
input parameters collected in Table 1. First of all, we would like to mention that we will 
not distinguish between the form factors of p and cu. Their difference is mainly caused 
by different values of the decay constants, /J^"*"^ 7^ fl^'^\ whose precise determination, e.g 
from experimental data for u — > e'^e~ , is complicated by mixing with the (f) meson. In the 
present paper we take the view that the uncertainty induced by letting T^^'' — Tf^'^ is 
negligible compared to other uncertainties. 

The longitudinal decay constants fpj^, can be extracted from the experimental decay 
rates r~ {p~,K*~)ut- as [22] 

fl = (0.209 ± 0.002) GeV, /J^. = (0.217 ± 0.005) GeV. 

There is no direct experimental measurement of the tensor decay constants f^^*, which 
instead have to be determined from non-perturbative methods such as lattice calculations 
[23, 24] or QCD sum rules [25, 18]. Lattice results are available in the quenched approx- 
imation with a chirally improved lattice Dirac operator, which allows one to reach small 
quark masses, and for the ratio of decay constants fy/fv [24]; a first study for the p with 
dynamical fermions was reported in Ref. [26]. One result of these calculations is that the 
ratio of decay constants only weakly depends on the quark masses. For the p, Ref. [24] 
quotes 




(2 GeV) = 0.72 ± 0.02 , 

latt 

obtained for the lattice spacing a = 0.15 fm. As for QCD sum rules, the value /^(l GeV) = 
(0.160 ± 0.010) GeV was obtained in Ref. [27]. For the present paper, we have re-analysed 
the corresponding sum rules, using updated values of ag and NLO evolution of f^, and 
find 

f^{l GeV) = (0.165 ± 0.009) GeV. (6) 

Also fp can be calculated from sum rules, yielding (0.206 ± 0.007) GeV. If one calculates 
the ratio directly from QCD sum rules, one finds^ 

^ I (2 GeV) = 0.69 ± 0.04 , 

■fp J SR 

in agreement with the lattice result. 

The determination of is less straightforward, see Ref. [18], where 

/^,(1 GeV) = (0.185 ± 0.010) GeV (7) 



iThe NLO scaling factor f^{2GeV)/ f^{l GeV) is 0.876. 
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NLO twist-3 
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■^3K* 


0±2 


TP 


NLO twist-3 


G-odd, UR 




AT 
S4p 
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[21] 
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LO twist-4 


UR 



Tabic 1: Hadronic parameters entering i?i,2,3 in the sum rule for Ti, Eq. (4). For twist- 
2 DAs, we use both the conformal expansion Eq. (8), truncated after n = 2, and the 
model of Ball and Talbot (BT) [19] given in terms of two parameters, A and p. The 
values of a2{p) and a2{K*) are highly correlated; we fix a2{K*) — a2{p) = —0.04 ± 0.02 for 
both longitudinal and transverse DAs. The twist-3 and -4 G-odd parameters have never 
been considered before; all twist-3 and -4 parameters are under revision (UR) and will be 
considered in full detail in Ref. [12]. In this paper (TP), twist-3 SU(3)-breaking effects are 
only taken into account at LO in the conformal expansion. Higher-orders in the conformal 
expansion of twist-4 DAs are calculated in the renormalon model of Ref. [17], see text. 
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was obtained. Evaluating the ratio fx*/ fx* directly from sum rules, we find 



(2GeV) = 0.73 ±0.04, 



K* / SR 



which agrees with the interpolation between the corresponding results for p and (f) obtained 
from lattice [24]. 

Summarising, it is probably fair to say that the present status of fy decay constants 
is not entirely satisfactory. The accuracy of the QCD sum rule estimates is unlikely to 
improve, so any significant reduction of uncertainty has to come from lattice. For the 
moment, however, all existing lattice results still come with considerable uncertainty (no 
continuum limit, no results for K* with chirally improved Dirac operator), so that in the 
numerical analysis of ^ we will use the experimental results for fj,^x* ^^'i QCD sum 
rule results (6) and (7) for f'^^*- 

As for twist-2 DAs, the standard approach is to parametrise them in terms of a few 
parameters which are the leading-order terms in the conformal expansion 

p^) = 6«(1 - + a„(/x2)C3/2(2« - 1) j . (8) 

To leading-logarithmic accuracy the (non-perturbative) Gegenbauer moments a„ renormal- 
ize multiplicatively. This feature is due to the conformal symmetry of massless QCD at 
one-loop, the start to mix at next-to-leading order, see appendix. Although (8) is not 
an expansion in any obvious small parameter, the contribution of terms with large n to 
physical amplitudes is suppressed by the fact that the Gegenbauer polynomials oscillate 
rapidly and hence are "washed out" upon integration over u with a "smooth" (i.e. not too 
singular) perturbative hard-scattering kernel. For vector mesons, one usually takes into 
account the terms with n = 1,2; the are estimated from QCD sum rules which are 
known to become less reliable for larger n. As an alternative, one can build models for 
based on an assumed fall-off behaviour of a„ for large n. The model of Ball and Talbot 
(BT) [19], for instance, assumes that, at a certain reference scale, e.g. p = 1 GeV, the a„ 
fall off as powers of n: 



oc 



{n + 1)P 

BT fix the absolute normalisation of the Gegenbauer moments by the first inverse moment: 



which can be viewed as a convolution with the singular hard-scattering kernel 1/u and 
gives all a„ the same (maximum) weight 1. The rationale of this model is that the DA is 
given in terms of only two parameters, p and A, and allows one to estimate the effect of 
higher order terms in the conformal expansion on observables. In this paper, we calculate 
the form factor using both conformal expansion, truncated after n = 2, and the BT model. 
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normalised by 02 and taking into account terms up to n = 8. We shall see below that the 
effect of terms with n > 2 is very small. 

For the p, a^'" have been determined in Ref. [27]. In the present study we have re- 
examined the corresponding sum rules and find, at the scale fi — 1 GeV, a^(p) = 0.15±0.07 
and a|(p) = 0.14 ± 0.06, which is shghtly smaller than the results quoted in Ref. [27]. As 
both values are nearly equal, we shall use a common value 

a^(p) = 0.15 ±0.07 = a|(p). (9) 

The corresponding value of A is 

= 1.24 ±0.11 = Aj, 

with a central value slightly larger than that used in Ref. [13]. The value of a2{K*) has been 
determined in Ref. [25]. Again, we re-examine these sum rules for the present paper. We 
find a2 (i^*) = 0.11 ± 0.09 and oiiK*) — 0.10 ± 0.08, which is more conveniently presented 
by the difference between 02 (i^*) and 02 (p): 

a^{K*) - a^{p) = -0.04 ± 0.02 , 

al{K*) - al(p) = -0.03 ± 0.02. 

As both differences are nearly equal, we shall use 

a2{K*) - a^ip) = -0.04 ± 0.03 (10) 

for both polarisations. This translates into A;^, = 1.18±0.14, with errors largely correlated 
with those of A^. 

The value of ai(p) vanishes by G-parity. The values of a\''^{K*) have been subject to 
some controversy over the recent years, which was settled only very recently; in this work, 
we use the values obtained in Ref. [18]: 

a^(X*) =0.04 ±0.03, a!(X*) = 0.03 ±0.02. (11) 

All odd Gegenbauer moments, i.e. the antisymmetric contribution to 0(m), can be re- 
summed using the same power-like behaviour of large moments as in the BT model. This 
model is also discussed in Ref. [19] and normalised to Oi; we include terms up to n = 9. 

Twist-3 and -4 DAs of vector mesons have been studied in Refs. [10, 11]. The results 
are complete for mesons with definite G-parity (with equal- mass quarks), but miss certain 
G-parity-breaking corrections. A complete analysis of all these corrections is in preparation 
[12]; here, we include those results that are already available [28]. In Ref. [10], the two- 
particle twist-3 DAs g^l'""^ and have been expressed in terms of integrals over the 
twist-2 DAs and the thrcc-particlc twist-3 DAs A, V, T. These integral relations are 
complete, but the explicit expressions for the three-particle twist-3 DAs given in [10] have 
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to be extended to include G-parity-breaking corrections as follows: 
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A{a) = SeOaictsa^Ca <j 1 + A^(ai - eta) + (^t { ^0,3 - - 



T(a) = 360q;iq;2Q;3 jreg + ^ C3 '^3 («i - "2) + «3 ^ -as - - 



(12) 



Here and ujs are G-parity conserving quantities, whereas ^3 and A3 are G-parity breaking. 
As K3 contributes to the form factor only at (!?(my/m^), and the A3 parameters are of 
non- leading conformal spin, we set, in the present analysis, the central values of all these 
parameters to zero and only take into account «;|. A QCD sum rule estimate of this 
parameter yields [12, 28] 

4(1 GeV) = 0.001 ±0.001. (13) 

The effect of non- zero values of k,^ and A3 is taken into account by the variation of these 
parameters around zero within the range given in Table 1; the dependence of on NLO 
G-parity breaking parameters is very small, as expected. 

The two-particle twist-4 DAs hs and A_l have been discussed in Ref. [11]; they are 
given by integrals over chiral-odd twist-4 three-particle DAs. The determination of the 
conformal-expansion coefficients of the latter is complicated by the fact that they contain 
"kinetic" mass-correction terms given by twist-2 matrix elements, which, to date, have not 
been obtained in a closed form, but have to be unravelled order by order in the conformal 
expansion. In addition, the direct determination of the "genuine" twist-4 corrections from 
QCD sum rules becomes increasingly complicated at higher-order conformal spin. For that 
reason, we invoke an alternative estimate of these corrections based on the renormalon- 
model developed in Ref. [17]. The general idea of this technique is to estimate matrix 
elements of "genuine" twist-4 operators by the quadratically divergent contributions that 
appear when the matrix elements are defined using a hard UV cut-off. In this way, three- 
particle twist-4 DAs can be expressed in terms of the leading- twist DA (f)± [17]: 



Ti{a) = -Tsia) 



CI 



;i-a2) 



(1 - a^Y 



(j)±{ai) 0_l(1 - «2) 



1 — «! 



1-^2 



T 



S{a) = -S{a) = -U 



0±(«i) _^ 0±(1 - a2) 
1 — cti 1 — q;2 



(14) 



The above formulas differ from those given in Ref. [17] by the change of argument a2 — > 
1 — q;2 in the second terms on the right-hand side; this is to properly account for G-parity- 
brcaking effects [29] . One prediction of the renormalon model is that the two independent 
LO twist-4 couplings and (J add up to 0, which is consistent with the direct calculation 
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Figure 1: Contribution of Ri to the convolution integral in (4) as a function of u. Solid 
line: long dashes: /|(mp/mft)i?2(M), short dashes: fj-{mp/mbfR^{u). 

from QCD sum rules [21]. The above formulas also allow one to estimate the "genuine" 
twist-4 G-parity breaking contributions to 7] and S; we refrain from giving explicit 
formulas in this paper, but refer to Ref. [12]. For the calculation of the contribution of 
twist-4 terms to ^, wc use two methods: firstly the full renormalon model (14), and the 
corresponding expression for Ax as given by the equations of motion [11] {hs = in this 
model). This accounts for the genuine twist-4 corrections; the "kinetic" corrections, as 
far as they are known, are added using truncated conformal expansion. Secondly, we use 
truncated expansion for all twist-4 DAs, describing G-parity-breaking terms by the values 
they assume in the renormalon model, see Refs. [29, 12] for more details. The predictions of 
both methods for the end-point behaviour of the DAs near u — 0,1 differ quite drastically; 
nonetheless, both prescriptions given nearly the same result after integration over u. 

One more parameter that enters the kinetic mass corrections to twist-3 and -4 DAs, 
induced by the equations of motion, are the quark masses „ ^- We choose rns{2 GeV) = 
(0.10 ± 0.02) GeV, which is in accordance with both lattice [30] and QCD sum rule calcu- 
lations [31], and let rriu^d = 0. 

With all DAs available, we can now assess the respective size of the contributions of 
the various Ri to the sum rule (4). To this end, we plot, in Fig. 1, the functions Ri for the 
p meson, multiplied by the corresponding weight factors, for u > 0.5 which is about the 
smallest value of uq. The plot clearly shows that Ri is dominant. It also shows that i?2,3 
exhibit (integrable) end-point singularities for — > 1. Based on these results, we expect 
the impact of the first neglected term in the light-cone expansion, which is 0(my/m^), to 
be very small. 

Before we can evaluate the sum rule for ^, we also have to discuss the choice of mj, and 
the sum-rule-specific parameters and Sq. The good news is that, although numerator 
and denominator of (5) both depend on rrih, Sq and M^, this dependence cancels to a large 
extent in the ratio. The reason hereof is quite evident from (4): controls the respective 
weights of contributions of different u; as these contributions are nearly equal in numerator 
and denominator of (5), except for moderately sized SU(3) breaking, it follows that one 
can choose equal in Tf^^^'^ ^ and that the resulting dependence on should be very 
small. This is borne out by the left panel of Fig. 2, where we plot ^ as function of for 
central values of the input parameters and Sq — 35GeV^. For comparison, we also show 
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Figure 2: Left panel: ,^ as a function of the Borel parameter for Sq = 35 GeV^ and 
central values of the input parameters. Right panel: ^ as a function of the continuum 
threshold sq for = 8 GeV^ and central values of the input parameters. Solid lines: DAs 
in conformal expansion; long dashes: BT model [19] for twist-2 DAs; short dashes: BT 
model for twist-2 DAs and renormalon model for twist-4 DAs [17] . 
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Figure 3: ,^ as a function of /^*(1 GeV). Solid line: /^(l GeV) = 0.165 GeV, dashed lines: 
shifted by ±0.009 GeV. 



1.3 

1.275 
1.25 

1 . 225 
1.2 
1 . 175 
1.15 
1 . 125 



0. 0.01 0.02 0.03 0.04 0.05 0.06 0.07 




0.08 0.1 0.120.140.160.18 0.2 0.22 
02 (p) 



Figure 4: Left panel: ^ as a function of a^iK*^ at 1 GeV. Right panel: ^ as a function of 
a2(p) at 1 GeV. Sohd hue: a2{K*) = a2{p) - 0.04; dashed lines: a2{K*) shifted by ±0.02. 
Longitudinal and transverse parameters a| and are set equal. 
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^ calculated using the BT model for the twist-2 DAs (long dashes) and, in addition, the 
renormalon-model for twist-4 DAs (short dashes). All three calculations agree with one 
another very well. The fact that the impact of the BT model is only minor shows that 
the sum rules are sensitive only to a few gross characteristics of the twist-2 DAs, but not 
to the details of their behaviour near the end-point u — 1. As for the renormalon-model 
DAs, the difference to the truncated conformal expansion is most marked for small values 
of M^, which can be easily understood from the fact that for smaller the weight of 
contributions from u close to 1 gets enhanced and hence the difference between the end- 
point behaviour of conformally expanded DAs and renormalon- modelled DAs becomes 
more visible. The right panel of Fig. 2 illustrates the effect of a variation of Sq for fixed 
M^, which is also small. The value of Sq sets the lower limit of the integral over u, and 
again the dependence on sq largely cancels in the ratio as the integrands arc equal up to 
SU(3)-breaking effects. As sq itself is nearly independent of the final-state meson, it is 
natural to choose the same value in both numerator and denominator. As for rrib, it only 
enters in the ratio mv/rrif,, which controls the respective contributions of -Ri,2,3, the lower 
limit of integration uq = ml/so and the Borel exponential exp(— m^/(iiM^)). In the latter 
two parameters, a change of mj, is effectively compensated by a change of sq or M'^, which, 
as we have just discussed, induces only very small variations of the sum-rule result. The 
ratio mv/rrib changes from 0.185 for K* and nib = 4.8 GcV to 0.193 for rrib = 4.6 GeV, 
which also has only very minor impact. Based on these observations, we choose to evaluate 
^ for fixed rrib — 4.8 GeV, = 8 GeV^ and Sq = 35 GeV^ and attach to ^ a corresponding 
uncertainty of ±0.005. 

We are now in a position to obtain a result for ^ and estimate its uncertainty. The 
dominant uncertainty is due to the dependence of ^ on the chiral-odd twist-2 parameters. 
In Fig. 3 we plot ^ as a function of /^*(lGeV), for various values of /^(IGeV). The 
uncertainty in both parameters causes an uncertainty in ^ of ±0.08. In Fig. 4, left panel, 
we show the dependence of ^ on ai{K*), which induces a change in ^ by ±0.03; the variation 
of a^{K*) and a\{K*) as separate quantities induces the same change. The right panel 
shows the dependence on a2 which is rather mild and causes ^ to change by ±0.02. The 
variation of the remaining parameters within the limits specified in Table 1 causes another 
±0.02 shift in ^, so that we arrive at the following result: 

C = ^ ^ = 1.17 ± 0.08(/^-^;^*) ± 0.03(ai) ± 0.02(02) ± 0.02(twist-3 and -4) 

(0) 

± 0.01 (sum- rule parameters, nib ^-nd twist-2 and -4 models) 

= 1.17 ±0.09. (15) 

The total uncertainty of ±0.09 is obtained by adding the individual terms in quadrature. 
Let us stress again that the error of this result is dominated by far by parameter uncertain- 
ties, and is nearly independent of the sum rule specific parameters; it is also independent 

Of/B. 
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3 Determination of |Vtd/Vts| 

Let us now turn to the calculation of the ratio of branching ratios and the determination 
of The Belle collaboration has measured the quantity 



R. 



exp — -p= 



where B{B {p,u)^) is defined as the CP-average ^[8(3 (p,uj)^) + B(B (p,uj)^)] 
of 

B{B (p, uj)^) = i ^ p+7) + ^ [B{B' - A) + B{B' - ^l)] } , 

and B{B — > K*'f) is the isospin- and CP-averaged branching ratio of the B 
channels. In 2005, Belle reported a 5.1a measurement [1], 



R^^^'"" = 0.032 ± 0.008(stat) ± 0.002(syst) , 
followed by a 5.2(T measurement by BaBar in 2006 [3]: 



R. 



BaBar 
exp 



0.024 ± 0.005 , 



(16) 
(17) 



where the statistical and systematical uncertainty are added in quadrature. HFAG com- 
bines both results into the average [32] 



i?HFAG ^ Q ^ 0.005 . 



(18) 



Within QCD factorisation, and using the notations of Ref . [6] , the amphtude for B 
can be written as 



A{B Vj) 



V2 



where Xg are products of CKM matrix elements and the factorisation coefficients a"'^ consist 
of Wilson coefficients and non-factorisable corrections from hard scattering and annihila- 
tion; explicit expressions can be found in Ref. [6]. ay'^ depends in particular on the form 
factor Ti and the twist-2 DA (j)v-±- The theoretical expression for R is then given by 



R 



th 



Vtd 




Vts 


e \ 



'pi "'B 



«7(P7) 



1 -|- Re {Sa± + Sao) 



+ i(|<5a±r + |5aon{ — 



2Ri, cos J + Rl 



Rl — Rb cos 7 
l-2RbCos-f + Rl_ 

(19) 



with 5ao,± = cijip'^'^l) /(^rip^'^l) — 1- Here, 7 is one angle of the CKM unitarity triangle 
and Rh one of its sides: 
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Equation (19) differs from the expression given in Ref. [6] by the terms in \5a\^ which were 
neglected in that paper. It is obtained in the SM, assuming that B{B^ p^'-f) = B{B^ 
a;7), which indeed should be the case up to a small difference in the decay constants, a 
tiny difference in phase space and the sign of the contribution of weak annihilation (WA) 
diagrams, which is also small numerically. Equation (19) is also valid in extensions of the 
SM where the CKM matrix is still unitary and the a/ do not carry a new weak phase, for 
instance Minimal Flavour Violation; in this case new physics could change the values of 
(5ao,±. 

Let us first discuss the dependence of (19) on CKM parameters, described by the terms 
in square and curly brackets. The up-to-date value of iKj^/Kfel, as provided by the heavy 
flavour averaging group HFAG in March 2006, is, adding errors in quadrature [32]: 



V, 



ub 



cb 



0.106 ±0.008. 



In Fig. 5 we plot the CKM factors 

— Rb two ) ± i,^ 



_ Rf-Rbcosj _ Rl 

2Rh COS 7 + 1 — 2Rh cos 7 + -R^ 



as functions of 7; the uncertainty induced by R^ is small. What is the currently preferred 
value of 7? HFAG is yet to provide averages of the individual results obtained by BaBar 
and Belle, so we use the value quoted by the UTflt collaboration in March 2006 [33]: 

7uTfit = (71±16)°, (20) 

which is obtained from tree processes only and hence can be assumed to be free of new 
physics. We then obtain 

/cKM = 0.07 ±0.12, ^CKM = 0.23 ±0.07. (21) 

As /cKM is rather small for the angle 7uTfit, Eq. (20), the contribution of the corresponding 
non-factorisable contributions to (19), collected in Re5ao,±, is heavily suppressed. 

The parameters 10^(^7) | and [07(^*7)1 are almost exactly equal, so we set \a'j{p'j)/ 
aj{K*^)\ — 1. Is that value hkely to be changed by non-factorisable corrections? One 
type of such corrections, soft-gluon emission from charm loops, has been calculated in 
Refs. [34, 28]; it amounts to a contribution to \aj\ of 0{1/ {rribml)) and is small by itself 
(~ 2%), and even smaller is its SU(3)-breaking that changes {ajlp'j) / aY{K*'~f) \ by less than 
1%. Another source of corrections comes from terms in (a"(-ft'*7) — a^{K*^)) / a^{K*'y) , 
multiplied by the CKM factor ~ A^/ckmj which is tiny indeed. Hence, we do not see any 
obvious source of significant corrections to \a^{p'j) / a^{K*'j) \ = 1. 

The value of (5ao,± is calculated in QCD factorisation, which is accurate to 0{as), but 
misses, in general, terms that are suppressed by inverse powers of m^^c- The most important 
of these power-suppressed corrections is weak annihilation, which can actually be calculated 
in QCD factorisation, at least at tree level. WA is CKM suppressed in B ^ K*'-/ and 
mainly affects 07(^^7), but it is small for a'^{{p^,u!)j) because of a suppression by Wilson 
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0. 25. 50. 75. 100. 125. 150. 175. o. 25. 50. 75. 100. 125. 150. 175. 

7 7 

Figure 5: Left panel: the CKM factor /ckm in Eq. (19) (square brackets) as a function 
of 7 for \Vub/Va,\ = 0.106 (solid line) and iKfe/Kbl = 0.106 ± 0.008 (dashed lines). Right 
panel: ditto for Qckm (curly brackets). 



coefficients and the fact that WA is proportional to the electric charge of the quark involved, 
namely the u quark for a1^{p^^) and the d quark for a1^{{p^, Although WA is formally 

power-suppressed, it gets enhanced in a1^{p^^) by large Wilson coefficients and the absence 
of 0{as) suppression that affects other non-factorisable corrections. Numerically, WA is 
actually as large as the leading (in l/m^,) non-factorisable terms. In view of the importance 
of this contribution, we treat WA in two different ways: firstly, by using the QCD-factorised 
expression given in Ref. [6]; and secondly, by using the results obtained from QCD sum 
rules on the light-cone [35, 36, 8]. 

The WA contribution to the amphtude of e.g. B" — > p~7 can be written in the following 
way: 

Gf , f ^ 1 



A{B- ^ p 7)wA = ^ A„ j^Ci + - C2 J {p''i\{du)v-A{uh)v-A\B ). 
In naive factorisation, the matrix element on the r.h.s. can be written as 

{p~'-i\{du)v-A{uh)v-A\B~) = {p~\{du)v-A\^) {l\{uh)v-A\B~) 

+ {p-^\{du)v-A\Q){Q\{vh)v-A\B-) . 

The second term on the r.h.s. has been shown to vanish in the chiral limit, sec Ref. [37] 
for more details, so we will focus on the first term. Corrections to naive factorisation are 
of 0{as)i which may also relax the chiral suppression of the second term. Neglecting the 
latter, we have, in the notations of Ref. [36] 

A\ 

{p-{pHq)\{du)v-A{ub)v-A\B-{pB)) = vO^^e^f) {F^e'^'^^'^eWpB,?, 

-iFA[e^^^\pB ■ q) - q^(e^'^ ■ Pb)]} ■ (22) 

^Equation (22) differs from the definition given in Ref. [36] by an overall sign. The reason is that in 
[36] the covariant derivative D^j^ = d^ — ieA^ was used, corresponding to a negative value of e. In order to 
keep Fv^A positive, we change the sign of the definition. 
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The form factors Fa^ can be calculated in QCD factorisation themselves; both are 
then equal, and to LO accuracy one has 

F^T^ ^ ^aT'' = %^ (23) 

with Qu = 2/3 the electric charge of the u quark and the first inverse moment of the 
S-meson DA (pB- 

Jo ? Ab 

Equation (23) agrees with the result obtained in Ref. [6] by direct calculation of the WA 
diagram. Corrections are either of 0{as), and have been calculated in Ref. [38], or they 
are suppressed by powers of l/rrii,. The dominant source of the latter comes from photon- 
emission from a soft u quark and has been calculated in Refs. [35, 36], together with the 
perturbative photon emission giving rise to (23). The emission of photons from a soft 
quark line is governed by the parameter x, the so-called magnetic susceptibility of the 
quark condensate, {Q\q(Tai3q\0) f — V 47rQ; QqXiQQ) Fap , which has been discussed in detail 
in Refs. [21, 39], together with higher- twist DAs of the photon. Its contribution is, in 
the heavy quark limit, suppressed by one power of l/rrib with respect to (23), but at 
finite quark mass its size is set by the dimensionless parameter x{'^'>^)/fB ~ 0.2, with 
x(lGeV) = (3.15 ± 0.3) GeV"^ [39], which is not really small. In calculating the WA 
contribution to 6a±, we will use both expressions for Fy^A- -^wA°^' ^1- (23), and F^^^^ 
from the QCD sum rule calculation, sec Ref. [8] for details. 

Let us first discuss 6ao, where WA is suppressed and can be neglected. Its dependence 
on hadronic parameters is controlled by the factor /s/(Tj^^''Ab); it also depends, to a lesser 
extent, on and the twist-2 DA (f)p-±. To estimate the uncertainty of Re5ao and |5aop, 
we set /b = (0.205 ± 0.025) GeV, which is an average of quenched and unquenched lattice 
calculations [40, 41] and QCD sum rule determinations [42]. We also use Tf = 0.27 ±0.03 
from light-cone sum rules,^ and Xb{1 GeV) = (0.46±0.11) GeV, obtained in Ref. [43]. This 
value supersedes the guesstimate = (0.35 ±0.15) GeV [44] used in previous calculations 
and agrees with the value (0.46 ± 0.16) GeV found in Ref. [45]. We evaluate all spectator- 
interaction contributions, that is those involving Ab, at the scale /x^ = m% — ml, which is 
of order ~ Aqcd^6 as advocated in Ref. [6], but by a factor 2 larger; this is motivated, 
in part, by the fact that the anomalous dimensions governing the renormalisation-group 
running of the Wilson coefficients are given for 5 flavours only in Ref. [4] and hence should 
not be used at scales as small as ( Aqcd"t.6) ^''^ ~ 1.5 GeV. We then need to evolve Xb from 
1 GeV to which can be done using the following evolution relation [46]: 

Ab'(/^) = Xb\^^o) {l + ^ (1 - 2aB(/.o))} , (24) 

where (TB(lGeV) = 1.4 ± 0.4 is given by an integral over the S-meson DA 0b and was 
estimated in Ref. [43] . We then have 

XbM = (0.51 ± 0.12) GeV. 



^This value, and in particular its error, is quoted from our previous paper in Ref. [13], but is in agreement 
with the evaluation of Eq. (4). 
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We can now cast most of the dependence of 6ao on hadronic input parameters into a 
dependence on XsifJ'h) only, varying it in the interval = (0.51+[!:i?) GeV. We also allow 
for 20% power-suppressed corrections to the leading (in l/rrib) non-factorisable corrections 
and find 

ReSao = 0.06 ± 0.02(Ab, /b, Ti) ± 0.06(C>(l/m6)) , 

\Sao\' = 0.014 ±0.004(Ab,/b,Ti)^°:°09((^(1M))- (25) 

Let us now turn to da^. Neglecting the effect of WA, one has 5a-j- = Sa^. Varying A^, 
Jb and Ti as before, and allowing for 20% power-suppressed corrections to leading non- 
factorisable contributions, we find 

Re(5a2^°^ = -0.19 ± 0.09(Ab, /s, Ti) ± 0.06(O(l/m6)) , 

ISaT^'^l' = 0.051°;°^(Ab,/b,T0^°:°?(O(1M)), (26) 
in QCD factorisation and 

Re5a2°°^^ = -0.06±0.04(SR) ±O.O6(0(l/m5)), 
|^^QCDSR|2 ^ Q_Q2 ± 0.01(SR)l°:°^(C»(l/m6)) , (27) 

using QCD sum rules for the WA contribution. The SR error reflects the dependence of 
the result on the QCD sum rule specific parameters and Sq and the value of x- 
Taking everything together, we have 

2 

[0.75 ± 0.11(0 ± 0.03(07'^ 7, Rb)] , 

2 

[0.75 ± 0.11(0 ± 0.02(07'", 7, Rb)] . (28) 

This result makes it clear that the theoretical uncertainty associated with 5a is small and 
that the error is dominated by that of ^ — the reduction of which is mostly a matter of 
more accurate (lattice and QCD sum rule) calculations, but is not affected by uncalculable 
l/rrib corrections. Within the quoted accuracy, the two different methods to calculate 
the WA contribution agree. We would like to stress here that it is precisely the CKM 
suppression of 5ao^± which also suppresses their uncertainties and renders the application 
of QCD factorisation to S — > viable. 

We are now in a position to obtain values for iVtd/Vtsj. Comparing (28) with the 
experimental results (16), (17) and (18), we get 

Belle 

= 0.207 ± 0.016(th) ± 0.027(exp) , 

BaBar 

= 0.179 ± 0.014(th) ± 0.020(exp) , 



dQCDF 
-^th 



V, 



td 



V, 



ts 



^QCDSR 



Vts 



Vts 



Vtd 

Vts 



16 



td 



Vt 



ts 



HFAG 



= 0.192 ± 0.014(th) ± 0.016(exp) . 



(29) 



These values can be compared with that following from 7 and the unitarity of the 
CKM matrix: 



Vt 



td 



Vt 



ts 



A(l + - 2i?6 cos 7)^/2 ^ Q ± 0.029 . 



(30) 



SM 



Both results agree well within errors. As (30) is obtained from tree-level processes only, it 
represents the "true" value of iV^td/Vt^l in the SM. 

A third determination of |Vid/Vts| can be obtained from B mixing. In the SM, we have 



^rud ruB^ PbBb^ 



Vt 



ts 



V 



td 



(31) 



The current world average for Am^ is Anid — (0.507 ± 0.005) ps ^ [22]. Am^ has recently 
been measured by the CDF collaboration [47], 



AtUs = 17.77 ± O.lO(stat) ± 0.07(syst)ps-^ , 



(32) 



with an accuracy that exceeds 5a significance. DO provided a two-sided bound at 90% CL 
[48]: 

17 ps~^ < Arris < 21 ps"^ . (33) 

The hadronic matrix elements in (31) are obtained from lattice simulations. The most 
up-to-date results for the decay constants have been obtained by the HPQCD group, using 
unquenched Uf — {2 + 1) configurations [49]: 

/b7/s. = 1-20(3)(1), 

where the first error is statistical and from chiral extrapolation and the second comes from 
"other uncertainties" [49]. The particular strength of this calculation is that light quark 
masses as small as m^/S could be reached, which implies that only a moderate extrapolation 
to the physical chiral limit is required. As for the ratio of Bb^^, the currently best result 
is obtained from unquenched n/ = 2 calculations (JLQCD collaboration [50]): 

BBjBB,^lM7{16){tf,), 

where the first error is statistical and the second systematic. In this calculation, the 
minimal light quark mass was nig = 0.5ms, which requires a more substantial extrapolation 
to the physical limit and is responsible for the large systematic uncertainty. A combination 
of both results yields [41]: 



f 

J^'Bs _ . 210 ('+47^ 
jBd^B^ 



(34) 
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where the errors have been added in quadrature. This procedure may be problematic as it 
combines resuhs with different systematic effects, but yields the most reliable unquenched 
result to date."^ From this, one finds 



Vts 



= 0.2060lH^^^(th) 



± 0.0007(exp) , 



(35) 



Am 



which is the result obtained by the CDF collaboration [47]. This value, too, agrees with 
the two previous determinations. Finally, one can compare our result also to the results 
of global fits of the unitarity triangle. The UTfit collaboration quotes, in September 2006, 
[33] 

Vtd 



Vt 



ts 



0.202 ± 0.008 , 



UTfit 



whereas CKMfitter gets 



Vu 

Vts 



0.201 



+0.008 
-0.007 



CKMfitter 



Again, all values agree within errors. 



4 Summary and Conclusions 

In this paper we have presented a new analysis of the form-factor ratio ^ = t^-^^* /T^~*'' 
from QCD sum rules on the light-cone, paying particular attention to the size of SU(3)- 
breaking effects. We have obtained 

i = 1.17 ±0.09; 

this value is nearly independent of QCD sum-rule-specific parameters and the error is 
dominated by that of the tensor decay constants f^j^,. A reduction of these errors by a 
factor of two would reduce the total uncertainty to ±0.06. The numerical values of these 
constants come mainly from QCD sum rules, partly from quenched lattice calculations. 
A determination from unquenched lattice calculations with reduced errors would be very 
desirable indeed. 

We then have analysed the non-factorisable corrections to R = B{B ^ {p,uj)^)/B{B ^ 
K*j) in the framework of QCD factorisation. The dominant power-suppressed correction 
comes from weak annihilation diagrams that mostly affect B"^ — > p^7. We have estimated 
these corrections both in QCD factorisation and using QCD sum rules, and find that the 
results agree within errors; we will present a more detailed discussion of power-suppressed 
corrections in a separate publication [8]. Our present best estimate of i?th is given in 
Eq. (28). We then extracted the ratio of CKM matrix elements iVid/Vtsls^v-y from i?exp 
obtained by BaBar and Belle, respectively, and averaged by HFAG, and find the values 

''A critical discussion of these lattice results, and their impact on the constraints on new physics from 
B mixing, can be found in Ref. [51]. 
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given in Eq. (29). Our results for this parameter agree well with all other determinations 
available from various sources as summarised in the previous section. They also agree with 
the value extracted from B mixing, using the new measurement of Am^ reported by the 
CDF collaboration. Presently, there is no indication for new physics to be inferred from 
these results. 
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Addendum to v3 

Please note that in the arXiv version v2 of this paper, which is identical with the published 
version JHEP 04 (2006) 046, we used the BaBar bound quoted in Ref. [2], R^^^^' < 
0.029 at 90% CL, which was combined, by HFAG, with the Belle measurement to R^^'^ = 
0.024 ±0.006 and resulted in \ Vtd/Vts\B^_i-y^ = 0.179 ± 0.014(th) ± 0.022(exp). These values 
have changed with the BaBar measurement oi B{B — > (p, a;)7) reported in Ref. [3]; the 
corresponding new result for |Vtd/Vi<j| is given in (29). 



Appendix: NLO Evolution of Twist-2 DAs 



To leading-logarithmic accuracy, the (non-perturbative) Gegenbauer moments in Eq. (8) 
renormalize multiplicatively as 



C(/^^) = ^^"'/^^*^an(/^^), 



(A.l) 



where L — as{n'^) / o-silJ^Q), — (33 — 2A^j)/3, and the anomalous dimensions are given 

by 



7f) = 8C^U(n + 2) + 



7e - 7 



4 2(n+ l)(n + 2)y ' 

7n^^°)-8C^(^^(n + 2)+7^;-^). 

To next-to-leading order accuracy, the scale dependence of the Gegenbauer moments is 
more complicated and reads 



^NLO/,,2\ ^ /',,2a pNLO I 



47r 



n-2 



5:a.(;.^)L-r/(2/.o)4i), (A.2) 



where 



k=0 



JO) , 



a^(/x^) - as(yWo) 
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with L = as{n)/as{iJ.o), Pi = 102 — (38/3)A^/; "jn^ are the diagonal two-loop anomalous 
dimensions, which have been calculated, for the vector current, in Ref. [54], and, for the 
tensor current, in Ref. [55]. The mixing coefficients d^^^, k < n — 2, are given in closed 
form in Ref. [53] ; these formulas are valid for arbitrary currents upon substitution of the 
corresponding one-loop anomalous dimension.^ 

For the lowest moments n = 0, 1, 2 one has, explicitly: 

7o ^, 7i 243 81 243 81 ' 

ui) 724 104 ^^ m) ±(i) 38044 904 ,^ , , 

7o^^ = -9- - ^ ^/ ' ^1 = 124 - 8Nf , 72^(^) = - ' (^-3) 

and 

11(1) _ 35 20 - 3/jo . .50/(9/3o)-n 
"'^^ ~¥50-9/3o^ ^' 

^ 28 16-3^ . _ ^40/(9/Jo)-i) . (A.4) 
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